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O ■ Abstract 

y^ • We develop methods to extract resonance widths from finite volume spectra of 

' , 1+1 dimensional quantum field theories. Our two methods are based on Liischer's 

qj I description of finite size corrections, and are dubbed the Breit-Wigner and the im- 

r^ • proved "mini-Hamiltonian" method, respectively. We establish a consistent frame- 

^ . work for the finite volume description of sufficiently narrow resonances that takes 

'k>( I into account the finite size corrections and mass shifts properly. Using predictions 

^ • from form factor perturbation theory, we test the two methods against finite size 

data from truncated conformal space approach, and find excellent agreement which 
confirms both the theoretical framework and the numerical validity of the methods. 
Although our investigation is carried out in 1+1 dimensions, the extension to physi- 
cal (3+1) space-time dimensions appears straightforward, given sufficiently accurate 
finite volume spectra. 

1 Introduction 

Two-dimensional field theories have attracted quite a lot of interest for decades, partly 
because they are considered to be theoretical laboratories for the development and testing 
of quantum field theory methods. The present paper can be considered as one more step 
along this line, taking up the issue of developing methods to extract resonance parameters 
from finite volume spectra. 

The framevi^ork of tv^o-dimensional field theory is particularly v^^ell adapted to this 
problem because non-integrable cases are increasingly better understood. In integrable 
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theories there can exist stable excitations whose decay into lower mass particles would 
in principle be allowed by their mass and conserved internal charges, but is prevented 
by the infinitely many integrals of motion underlying integrability. However, by adding a 
perturbation that breaks integrability one can obtain a much richer phenomenology. At the 
same time, the underlying integrable model provides us with very powerful tools: firstly, 
the exact spectrum, scattering theory and form factors at the integrable point are known, 
and secondly, form factor perturbation theory (FFPT) [1, 2] can give a prediction for the 
mass shifts, phase shift corrections and decay widths when the non-integrable perturbation 
is switched on. Therefore the phenomenology is under firm control, which makes these 
theories into a very convenient testing ground^ 

On the other hand, there exists a rather remarkable tool to compute finite size spectra 
in two space-time dimensions, known as truncated conformal space approach (TCSA), 
originally developed in [4] and carried further in many papers, in particular [5, 6, 7]. 

Recently, Delfino et al. [2] gave a detailed description of the FFPT framework for the 
decay widths and performed calculations for the Ising model. They also proposed a way 
to extract the decay width from the finite size spectrum. However, the theoretical relation 
to the finite volume spectrum was determined incorrectly as we show in section 3.2.1. In 
addition, they could not reach the required precision to extract any quantitatively useful 
information from the TCSA spectrum either. In the course of our investigation we also 
present an improved realization of their proposal which is consistent with a Breit-Wigner 
resonance description, and demonstrate its quantitative validity using a more advanced 
TCSA analysis. In addition, we extend our considerations to another prominent example 
of non-integrable field theory, the double (two-frequency) sine-Gordon model. 

We set out to formulate and test methods to extract (partial) decay widths of unstable 
particles (resonances) from finite volume spectra, using the ideas in [8] as a starting point. 
In particular we develop new and efficient methods to extract the tiny effects of a narrow 
resonance with high precision, and test them in the framework provided by two-dimensional 
quantum field theories. We examine the scaling Ising model with a thermal and magnetic 
perturbation (i.e. a continuum version of the familiar two-dimensional Ising lattice away 
from critical temperature and in a magnetic field) and the double (two-frequency) sine- 
Gordon model, since these are the most studied examples of non-integrable quantum field 
theories in two dimensions (see the references given in sections 2.2 and 2.3). However, 
the methods we develop to extract the decay widths are very general, and their higher- 
dimensional extension seems to be rather straightforward. 

The paper is organized as follows. Section 2 introduces the framework of form factor 
perturbation theory and gives the theoretical predictions for the two models considered. 
Section 3 is devoted to the theoretical development of the methods used to extract the decay 
widths from finite volume spectra, and examines the consistency between them. Section 

^It is also possible for resonances to occur in an integrable field theory, as e.g. in the homogeneous 
sine-Gordon models [3], but the finite volume spectra of these models are not sufficiently understood for 
the purposes of the present investigation. In addition, the framework of perturbed integrable field theories, 
in which the decay can be switched on with the integrability breaking perturbation, provides us with a 
very useful control parameter, i.e. the perturbing coupling, which can be used to tune the resonance width. 



4 is devoted to the specification of numerical methods and extrapolation procedures used, 
and contains a discussion of the sources of numerical errors. We present and discuss our 
numerical results in section 5, and conclude briefiy in section 6. Two technical details are 
relegated to separate appendices: one is the FFPT derivation of the relevant decay width 
in the double sine-Gordon model and the other is the detailed calculation performed to 
examine the consistency of the methods developed in section 3. 

2 Unstable particles and form factor perturbation the- 
ory 

2.1 General formalism 

In the framework of conformal perturbation theory, we consider a model with the action 

A{fj., A) = .4cFT - At / dtdx^{t, x) - \ dtdx'^{t, x) (2.1) 

such that in the absence of the coupling A, the model defined by the action A{fJ,, A = 0) is 
integrable. The two perturbing fields are taken as scaling fields of the ultraviolet limiting 
conformal field theory, with left/right conformal weights /i$ = /i$ < 1 and hq, = h\^ < 1, 
i.e. they are relevant and have zero conformal spin, resulting in a Lorentz-invariant field 
theory. 

The integrable limit A{fi, A = 0) is supposed to define a massive spectrum, with the 
scale set by the dimensionful coupling /x. The exact spectrum in this case consists of some 
massive particles, forming a factorized scattering theory with known S matrix amplitudes, 
and characterized by a mass scale M (which we take as the mass of the fundamental particle 
generating the bootstrap), which is related to the coupling /x via the mass gap relation 

H = kM^-^^-^ 

where k is a (non-perturbative) dimensionless constant. 

Switching on a second independent coupling A in general spoils integrability, deforms 
the mass spectrum and the S matrix, and in particular allows decay of the particles which 
are stable at the integrable point. One way to approach the dynamics of the model is form 
factor perturbation theory initiated in [1]. Let us denote the asymptotic states of the A = 
theory by 

\A,i'^^)...A„{l^n))x=o 

which describe in states if the rapidities are ordered as t^i > . . . > ^„ and out states for 
^1 < . . . < t9„. Then the essential input to form factor perturbation theory consists of the 
matrix elements of the local field \E' (so-called form factors): 

Fl^,^ {A, ...,^n) = m{0, 0)|A,, (^i) ...A„ (^„)),=o 



which (due to integrability at A = 0) can be obtained in a closed exact form solving the form 
factor bootstrap axioms and identifying the appropriate solution using known properties 
of the operator \l/ (for a review of the form factor bootstrap, see [9] and references therein). 
Given these data the following quantities can be calculated to first order in A: 

1. The vacuum energy density is shifted by an amount 

5£,ac = A (0| ^ |0);,^o . (2.2) 

2. The mass (squared) matrix M'^^^ gets a correction 

5Ml, = 2\Fl{m,Q)5rn^,^, (2.3) 

(where the bar denotes the antiparticle) supposing that the original mass matrix was 



diagonal and of the form M'^^^ = ml6, 



ab 



3. The scattering amplitude for the four particle process a + b — > c + d is modified by 

rriamh smh v 

It is very important to keep in mind that this gives the variation of the scatter- 
ing phase when the center-off-mass energy (or, the Mandelstam variable s) is kept 
fixed [1]. Therefore, in terms of rapidity variables, this variation corresponds to the 
following: 



where 



d^ dX 

m.,f)m,. -\- m.. 

5^ = 



A=0 

maSnia + maSnia + {mhSnia + maSmh) cosh -(9 



niami, sinh § 
is the shift of the rapidity variable induced by the mass corrections (2.3). 

It is not yet clear how to extend these results to second order. However, it is possible 
to calculate the (partial) decay width of particles [2]. Suppose that the decay of particle 
Ac into particles Aa and Ai, is kinematically allowed: rric > ma + rrib. Then in the rest 
frame of particle A^. the rapidities of the outgoing particles are uniquely determined from 
energy-momentum conservation 

m„ sinh ^i""^)+m6 sinh ^['"^^ = 

rUa cosh i!)^^"-^' + mh cosh 'dj^"' = rric 

and the partial decay width can be calculated as 



;s^22l-<5„i 



cab [ ^^1 '^a 5 '^b 
c~*ab 



mllTLa 



sinh §a 



(cab) 



(2.5) 



All the masses nia in the above formulae correspond to those in the unperturbed (A = 0) 
theory. 
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2.2 Thermal perturbation of Ising model with magnetic field 

The action 

-Aising ih,T) = A^^i — h dtdxda{t,x) — t dtdxe{t,x) (2.6) 

describes the two-dimensional scaling Ising model, where A^^i is the action of the c = 1/2 
conformal field theory (free massless Majorana fermion), describing the conformal limit 
of the two-dimensional Ising model, e is the energy density operator (which is a primary 
field of weight A^ = A^ = 1/2) and a is the magnetization operator {Aa- = A^j = 1/16). 
The coupling h corresponds to the applied external magnetic field and r parameterizes the 
deviation of the temperature from criticality^. This model is considered as a prototype non- 
integrable field theory in [1] and is widely discussed in the literature [11, 12, 13, 14]. There 
are two options to consider it as a perturbed integrable field theory. One way is to take h 
as the perturbing coupling: in this case the Majorana fermions of the thermally perturbed 
Ising model can be demonstrated to go through some sort of confinement (McCoy- Wu 
scenario [15]), and the decay widths of the high-lying "meson" states were computed to 
0{h?) in [16]. Here we take the opposite route and consider r as the coupling that breaks 
integrability. 

For r = the spectrum and the exact S matrix is described by the famous E^ factorized 
scattering theory [17], which contains eight particles Aj, i = 1, . . . , 8 with mass ratios given 
by 



777,2 = 


TT 

= 27r7i cos — 
5 


7773 = 


TT 

= 27r7i cos — 
30 


7r74 = 


9 ^^ 

= 27^70 COS 

30 


7775 = 


277 
= 27^70 COS 

15 


7776 = 


TT 
= 27770 COS 

30 


7777 = 


= 27774 cos — 
5 


7778 = 


vr 

= 27775 cos — 

5 



and the mass gap relation is [18] 

mi = (4.40490857... )|/i|^/^^ 



or 



h = KhTu^^^ , Kh = 0.06203236 . . . (2.7) 



^Note that for this model the notation for the couplings is changed from ^, X to h,T in order to conform 
with the notations used in [2]. In addition, $ is now a and \l/ is e. 



We also quote the scattering phase shift of two Ai particles for A = 0, which has the form 

sinh i) + i sin ttx 



Sni^) 



1 
15 



{x} 



sinh {} — i sin nx 



(2.8) 



All the other amplitudes Sab are determined by the S matrix bootstrap [17]. 

The non-integrable model (2.6) has a single dimensionless parameter, which we choose 

as 

t 



T 



|/^|8/15 

The form factors of the operator e in the Eg model were first calculated in [19] and their 
determination was carried further in [2]. Using the exact form factor solutions the following 
results were obtained in [2] for the mass shift of Aa and the partial decay width associated 

to A^^Aa + Ab-. 
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= (-17.8933 


■ ■)(e) 
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■ ■)(e) 
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■ ■)(e) 


/4II 


= (36.73044. 
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/5II 


= (19.16275. 
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and the vacuum expectation value of the perturbing field e for r = is given by [20] 

{e) = eh\h\'/'' , e, = 2.00314... 
The dimensionless decay width for the processes Ac ^ Ai + Ai, c = 4, 5 can be written as 



mi 



rric 
nil 



4,5 



(2.9) 



2mi 



-1 



where 



/411 
/511 



36.73044 , 
19.16275, 



2.3 Double (two- frequency) sine- Gordon model 

Double sine-Gordon theory is another prototype of non-integrable field theories which 
can be understood by application of techniques developed in the context of integrable 
models [10]. It has several interesting applications e.g. to the study of massive Schwinger 
model (two-dimensional quantum electrodynamics), and in the description of a generalized 
Ashkin-Teller model (a quantum spin system), both of which are discussed in [10]. Another 
application to the one-dimensional Hubbard model is examined in [21] (together with 
the generalized Ashkin-Teller model mentioned above). A further potentially interesting 
application of the two- (and multi-) frequency sine-Gordon model is the description of ultra- 
short optical pulses propagating in resonant degenerate medium [22]. Its phase diagram 
was studied using non-perturbative finite size techniques in [23], and the results were 
recently extended to the multi-frequency generalization in [24]. There has been a certain 
doubt concerning the validity of factor perturbation theory (or rather the proper way 
of performing it) after Mussardo et al. [25] applied a semiclassical soliton form factor 
technique developed by Goldstone and Jackiw [26], and obtained results that contradict 
explicitely some of the results derived from standard form factor perturbation theory in 
[23]. However, the conclusions drawn from the semiclassical technique were shown to be 
untenable by extensive numerical work in [27], which upheld the results of form factor 
perturbation theory as applied in [23], and therefore we follow the same approach in this 
paper. 

The action of the two-frequency sine-Gordon model (with frequency ratio 1 : 2) is 

A = Ac=i + /U / dtdx cos (3(p + A / dtdx cos I —(p + 5 j (2.10) 

f 1 2 

Ac=i = / dtdx - (dip) 

The spectrum at A = consists of a soliton doublet with mass M related to /x via [28] 

/i = «:sg(OM^/«+^) 

2/(5+1) 




and breathers i?„ with masses 

m„ = 2Msin^ , n = 1, . . . , [1/^] (2.12) 

The exact S matrix of all these particles was derived from the axioms of factorized scat- 
tering in [29]. 

The strength of the non-integrable perturbation can be characterized using the dimen- 

sionless ratio 

_ A 

'' ~ 4+3; 

M2+2« 



where M is the soliton mass at the integrable point A = 0. From (2.12) 

— = l + 2cos7r^> 2 if ^ < 1/3 
nil 

so whenever B^ exists, its decay to 2 Bi-s is always kinematically allowed (for A small 
enough) and the same holds for all the -B„, n > 3. However, since 

— = 2 cos — < 2 

mi 2 

B2 cannot decay to a pair of -Bi-s (and therefore cannot decay at all) if the perturbation 
A and correspondingly the mass shifts are small enough. The decay width corresponding 
to the simplest process B-^ ^> Bi + Bi is calculated in Appendix A for 5 = —n/2. The 
dimensionless decay width can be written as 

^''^^^ - e ^ '^^ (2.13) 






I M ) 



M y \ 2mi 



where S311 is given by (A. 7). Because the B^ state is the lowest one whose decay can 
be observed, it is preferable to use in the comparison with truncated conformal space to 
minimize the truncation errors. In addition, there is a Z2 symmetry at 5 = — 7r/2 which 
can be turned to a significant numerical advantage (see the analysis in subsection 4.1 for 
details). 

3 Resonance widths from finite volume spectra 

In this section we discuss two different approaches to linking the finite volume spectra with 
the theoretical predictions of form factor perturbation theory in the previous section. The 
first one uses a Breit-Wigner parameterization for the finite volume levels in the vicinity 
of the volume Lq where the one particle level A^ crosses the two-particle level AaA^ before 
switching on the coupling. This method relies on obtaining the two particle phase shift 
from the finite size spectrum using Liischer's results [30] and then finding the Breit-Wigner 
resonance parameters directly from the phase shift. The second approach uses a "mini- 
Hamiltonian" description for the levels to extract directly the form factor determining the 
decay width. The two approaches are then compared and we learn how to improve the 
results of Delfino et al. [2] for the "mini-Hamiltonian" to make them consistent with the 
Breit-Wigner description (for a narrow resonance). 

3.1 Liischer's formulation: methods to extract the decay rate 
from the finite size spectrum 

Suppose that the finite volume spectrum of the theory is given in a form of functions Ei{L) 
describing the energy of the iih excited stated as a function of volume L, normalized to 
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the ground state energy in the same volume. Suppose further that we investigate a decay 
process Ac —* Aa + Af, {rric > rria + rrib) and that we have a parameter A describing the 
interaction responsible for the decay, such that for A = the particle Ac is stable (the 
notation here follows the conventions of section 2.1). 

We can then find the energy level corresponding asymptotically to a stationary particle 
Ac at A = 0, denoted by Ec{L) and two-particle levels composed of Aa and Ab with 
zero total momentum, denoted by Eab{L) (where we suppressed an index labeling states 
corresponding to different relative momenta of a and h). According to [30], up to corrections 
vanishing exponentially with L (which we neglect from now on), the finite size correction 
to the two-particle level is determined by the equations 

niaL sinh da + 5ab i'&a - A) = 2na7r 

nibL sinh ^^ + 6ab {^b - ^a) = 2n{,7r (3.1) 

and 

Eab{L) = iTLa cosh da + nib cosh db 

where 5 ah {d) = Sab {S) is the elastic two-particle phase shift defined from the elastic 
two-particle scattering phase Sab via the relation 

Sabm=e''^''^'>'^ 

and n = Ua = —rib for zero total momentum labels the different possible values for the 
relative momentum. As a result, the energy of any two particle level decreases as l/L"^, 
while the finite size corrections of a one-particle level are known to vanish exponentially 
fast with increasing volume [31]. In the infinite volume limit Eab{L) < Ec{L) (since in the 
limit L ^ oo they tend to m^ + nib and rric, respectively) and therefore there is a volume 
Lq where the two levels meet: Eab{Lo) = Ec{Lq) = Eq. 

3.1.1 The "mini-Hamiltonian" 

Switching on a small A 7^ the levels only move a small amount and thus can be easily 
re-identified, however the degeneracy at Lq is lifted by an amount related to the width of 
the decay process (as we shall see shortly). Let us restrict ourselves to the two-dimensional 
subspace spanned by the two levels. We suppose that the relevant region in the volume 
L is of order A, which turns out to be a self-consistent assumption. Then for L ^ Lq and 
A ~ the effective "mini-Hamiltonian" in this subspace which determines the behaviour of 
the levels to leading order m L — Lq and A can be parameterized as follows 

where the relative phases of the two states were chosen such that B is real and nonnegative 
(A, C are real by hermiticity). According to the prescription of degenerate perturbation 
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Figure 3.1: Behaviour of levels at A = and A 7^ 0, illustrated using actual numerical data 
for the Ising case. (Energy and distance are measured in units given by the value that the 
mass of the lightest particle (mi) takes at A = 0.) 

theory, the splitting of the two levels is given by diagonalizing (3.2) in the subspace of the 
two levels: 



SE{L) 



((ai - a2){L - Lo) + \{A{L) - C{L))f + AB^X' 



Neglecting the volume dependence of A, B and C for the moment, the minimum splitting 
between the two levels occurs at 



L/„ 



Ln 



-A 



A-C 
ai — 0:2 



(3.3) 



which shows that it was consistent to suppose that the relevant values oi L — Lq are of the 
order A. The minimal energy split is then 



6E{L, 



2\BX\ 



(3.4) 



This shows that to first order in A it is consistent to replace the coefficients A{L), B{L), C{L) 
by their values taken at Lq. This observation will be used several times in the sequel. 
The behaviour of levels discussed above is illustrated in figure 3.1. 

3.1.2 Breit-Wigner analysis 

Here we build on the approach advocated by Liischer in [8], but with some modification to 
make it more suited to numerical analysis. At A 7^ there are no stable one-particle states 
of Ac anymore, and therefore asymptotically all the states are identified as two-particle 
states AaAb. Due to the lifting of the degeneracies in the crossings with the one-particle 
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The AgAfj levels at X*0 
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(b) A 7^ 



Figure 3.2: Illustrating the behaviour of two-particle states in the integrable case A = 
and with a resonance present (A 7^ 0). 
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level of Ac these states acquire plateaux as illustrated in figure 3.2, which correspond to 
the appearance of a Breit-Wigner resonance contribution in the two-particle phase shift: 

6ab{E)=6o{E) + 6Bw{E) 

(as a function of the center-of-mass energy E), where Sq is the background part of the 
phase shift, and the Breit-Wigner contribution has the form 

, , E-Ec-iT/2 

with Ec as the center and F as the width of the resonance. If 

d6o 



^ iE 



<0 

E=Ec 



then the total phase shift has two extrema (a minimum oX E < Ec and a maximum at 
E > Ec, see figure 3.3). Their positions can be calculated to first order in V: 



Ei 
with the values 




6ab {E±) = 7c + 6o{Ec) ± (tt - 2 V'^ 



Now let us take two neighbouring two-particle levels Ei{L)and E2{L). Using the description 
(3.1) we can define the phase shift functions for the two levels as follows: first solve 

E = A/p^~+mf + Wp2 -\- ml -^ p{E) 
and then extract the phase shift as 

6,{E) = -Li{E)p{E) 

S^iE) = -L2{E)p{E) (3.5) 

where Li2(-E) are the volume-energy functions for the two states (here we supposed that 
the two-particle state has zero total momentum, and the momenta of the particles are p 
and —p). The functions 5{E) are illustrated in figure 3.3. Since 

5^{E) = 5ai,iE)-2miT 
HE) = 6ab{E) -2n27r 

where rii and n2 are the momentum quantum numbers which for neighbouring states satisfy 
n2 = ni + 1, we obtain 



A5 = mm5i{E)~max52iE)=A^/^ (3.6) 
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Figure 3.3: Phase shift functions Sn{E) extracted from the two-particle levels illustrated 
in figure 3.2 around the resonant energy Ef. . 

This gives a method to extract the value of T from the finite volume spectrum: determine 
first the phase shift functions Si^2 from two neighbouring two-particle levels and then take 
the difference between their extrema to find vT. Since T is of order A^, it follows that vT 
is of order A which means that it is much easier to measure than T in the small A regime 
when the resonance is narrow. This is a huge advantage over previously proposed methods 
which gave T directly, like the approach advocated in [8] which relates T to the slope of 
the plateaux at their middle point (more precisely at the value of the volume when the 
resonant contribution to the phase shift S of the given level passes through n) . Indeed it 
is already apparent from figure 3.2 that the level splittings determining A are much easier 
to observe compared to the slope at the middle of the plateaux. In addition, A can be 
extracted by measuring the spectrum in a small neighbourhood of Lmjn which makes the 
residual finite size effects (those that decay exponentially on the volume) much easier to 
control. 

Finally we remark that while the simplest possibility of extracting the resonance pa- 
rameters would seem to be the method of fitting the numerically determined phase shift 
function S{E) with a Breit-Wigner resonance function, in practice this has the severe dis- 
advantage that the finite size data necessary to extract the relevant part of the phase shift 
cover an extended range between two line crossings together with a neighbourhood of the 
crossings themselves. Over such an extended range, the contribution of residual finite size 
corrections varies substantially, which causes such a significant distortion of the shape of 
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Figure 3.4: Illustration of the residual finite size eff'ects in the Ising model, for the levels 
pertaining to A4 and two-particle states AiAi. The volume dependence of the two-particle 
levels follows very precisely the description given in eq. (3.1), but the variation of A4 
is entirely due to the residual finite size effects that we neglected. The plot shows the 
integrable case t = when the line crossings are exact, but the only result of switching on 
a small value for t is a slight shift in the lines and resolution of the degeneracy at the level 
crossings. 

the Breit-Wigner resonance contribution that there is no way to perform a reliable fit to 
the data. This problem is illustrated in figure 3.4, where it is obvious how much the one- 
particle level varies between two adjacent level crossings. This also makes the application 
of Liischer's "plateau slope" method impossible, which can be seen from the fact that in the 
absence of residual finite size effects we would expect the slope to be negative, as illustrated 
in figure 3.2. However, it is obvious from the data in figure 3.4 that the leading residual 
finite size corrections are such that the particle masses approach their infinite volume from 
below, which is known to be a generic feature in two dimensional field theories [33]. 

3.1.3 Linking the "mini-Hamiltonian" with the Breit-Wigner formula 

We can now match this method with the previous one. For simplicity we shall assume 
that the two decay products are identical to a particle denoted by Ai (i.e. a = b = 1). 
Parameterizing the two levels using (3.2) we can express the parameters in the "mini- 
Hamiltonian" by extracting 

A6 = rain6i{E) — m.ax62{E) 

from the levels and matching it with the result (3.6) of the Breit-Wigner analysis. From 
(3.2), the two levels read 

^1,2 = ^0 + ^^Y^^ + ^ («i + «2) {L - Lo) ± ^y/{ai-a2y{L-L^,,,y + 4B^X^ 



14 



1,2 - 


= mc + AX- 


-a{L- L,„in) ± a/«2(L - 

r T A^-^ 


-L„,„)2 + 52A2 


(3.7) 
(3.8) 



(E1/E2 corresponds to the sign choice +/— , respectively) where Lmin is given by eq. (3.3). 

At A = 

1 1 

£^1,2 = Eo + -{ai + 02) {L - Lo) ± - (ai - ^2) 1-^ - -^ol 

Recalling that nic is the mass of the particle A^ whose decay process we are interested in, 
we get Eq = rric (where we neglected the residual finite size eff'ects decaying exponentially 
with the volume, as before). Furthermore, it is known that E2 = rric for L < Lq and 
El = rric at L > Lq (due to the level crossing at Lq, the interpretation of the levels in 
terms of asymptotic particle states swaps over). As a result, ai = and 02 = —2a < , 
and so the expression for the two levels turns into 



with 



and denoting 

A = AiLo) , B = BiLo) , C = C(Lo) 

since as pointed out in section 3.1.1, the coefficients A{L) , B (L) , C (L) can be substituted 
with their values taken at L = Lq if one calculates only to the lowest order in A. 

Note that for \a{L — Lmm)! ^ \B\\ one of the levels (the identity of which depends on 
the sign of L — L^rim) is ^-independent and takes the value 

Ec (A) = rric + AX (3.9) 

which therefore can be identified with the mass of the resonance corresponding to the 
unstable particle Ac, to first order in A (this shift is also illustrated in figure 3.1). 
Using form factor perturbation theory (2.3) we obtain the identification 

A=— , fcc = Fcc{r7T,0). 

rric 

After some calculation (the details are given in Appendix B) the phase shift difference 
can also be extracted 



min(5i — max(52 = 2B\^I y^ — 4mf —— 



E=Eo 



where (5i 2 are the phase shift functions for the levels Ei^2{L) as defined in (3.5). Now this 
must be compared to the result from the Breit-Wigner analysis (cf. eq. (3.6) ) 



min 61 — max 82 




E=Eo 
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This way we obtain the relation 



2\2 



4a 



yml — 4ml 



On the other hand, using the form factor perturbation theory result (2.5) gives 

x2 
p _ -12 Jell 



mimi 



sinh di 



(cU) 



and so 



B 
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fell — 


F!u{^rr,4''\^r') 




rric = 2mi cosh 1!)^^ 


fell 


1 8a 



rrir 
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We also need an expression for a. At A = we can express it using the energy Eii{L) of 
the AiAi two-particle state from (B.2) 



IdE 



a = — 
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2 dL 



L=Lo 



Writing the energy of the two-particle level with the one-particle momentum p as Eii{L) 
2iyp(Lp + mf and using Eii{Lo) = m,c + 0(A) we get 



B 
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cll 



m, 
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ldp{L)\ 
p dL J 



(3.10) 



L = Ln 



Putting B = the level corresponding to the two-particle state AiAi takes the following 
value at Lq: 

En (Lo) = me + CA 

so C\ is the energy shift of the two-particle level at the volume L = Lq if the resonance 
contribution is neglected. This is also consistent with eq. (3.8) which for _B = gives the 
location of the level crossing to first order in A. This observation provides a possibility to 
determine C (for details see Appendix C.l), with the result 



C 



ldp{L)\ 
p dL J 



L = Ln 
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(3.11) 



3.2 Theoretical determination of the "mini-Hamiltonian" 

3.2.1 Mismatch between the Breit-Wigner formalism and the DGM predic- 
tion 

Now let us turn to theoretical determination of the "mini-Hamiltonian" (3.2) 
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H = Eo + {L- Lo) 



«! 



0:2 



+ A 



A{L) B{L) 
B{L) C{L) 



It was shown in subsection 3.1.3 that matching the "mini-Hamiltonian description" with 
the Breit-Wigner analysis gives the relations 



fee 
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rrir 
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ldp{L)\ 
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(3.12) 
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In addition, matching the Liischer description (3.1) of the two-particle states with the 
"mini-Hamiltonian" also gave (3.11) for C. 

Delfino, Grinza and Mussardo (DGM for short) obtain the following results for the 
matrix elements [2] 
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(3.13) 



where 

/,, = F| (zvr , 0) 

and 
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(3.14) 
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(3.15) 



(recall that the phase of the matrix element B{V) can be transformed away by redefining 
the relative phase of the two states, so we can take the absolute value in (3.15)). 
We can see that the result for A{L^ matches perfectly, but B{Lq) agrees only if 



\_dp{V)\ 
p dL J 



L=Lo 



1 
To 



This is only true if we neglect the background phase shift ^o since the quantization rule for 
the momentum is 



Putting 6q = 



p{L) 



pL + 5o {E) 


= 2mT 






E 


= 2^Jp^ + ml 




2mi 

L ^ 


( ^dp{L)\ 
\ p dL ) 


L=Lo 


1 



However, the background phase shift is already present at A = and therefore the formulae 
(3.13) are not consistent with the Breit-Wigner analysis to first order in A. In fact, the 
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background phase shift (especially the fact that Sq{Eq) < 0) played an important role in 
the considerations based on the Breit-Wigner resonance parameterization. 

In addition, the DGM prediction for C does not agree with (3.11) and so is inconsistent 
with the behaviour of two-particle states in finite volume even in the absence of a resonance. 
Taken together, this means that the DGM formulae (3.13) cannot be used self-consistently 
to extract the form factors from the behaviour of the levels for small A. 

3.2.2 Rederiving the "mini-Hamiltonian" from the Lagrangian approach 

We shall now resolve the above problem by reexamining the derivation of the "mini- 
Hamiltonian" from the Lagrangian approach. We normalize the one-particle states in 
infinite volume as 

{Aa{Pa)\Ab{pb)) = 2'JT6abS{Pa - Pb)Ea 

In finite volume the one-particle momenta (up to corrections vanishing exponentially with 

the volume) are quantized as 

27rna 

Pa=^^ 

with the density of states 

dn L 

dp 27r 
therefore the normalization reads 

{Aa{Pa)\Ab{pb)) = 5ab5paPbLEa 

To obtain an orthonormal basis the one-particle states in finite volume must be related to 
the ones in infinite volume as 

\Aa(j>a))L = -y==\Aa{pa)) 

The normalization of two-particle states in infinite volume reads 

{Aa{Pa)Ab{pb)\Ac{Pc)Ad{pd)) = 2n5ac5{Pa - Pc)Ea2TlSbdS{pb - Pd)Eb 

Taking the density of the two-particle states to be the product of the one-particle densities, 
the following expression is obtained for the states normalized in finite volume L 

\Aa{Pa)Ab{pb))L = -y==^==\Aa{pa)Ab{pb)) 

In particular, considering states containing two identical particles and having zero total 
momentum we get 

\Aa{Pa)Aa{-Pa))L = "^ | Aa(Pa) Aa(-Pa)) 
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where E = 2Ea is the total energy of the state. The perturbing operator reads 



Jo 



The mini-Hamiltonian can be computed taking the matrix elements of H' in the basis 

|^c(0))l=Lo , \Aiip)Ai{-p))L=Lo 

Using translational invariance to perform the spatial integration, and the fact that the 
energy of these states is Ec{Lq) = Eii{Lq) = rric the expression of the "mini-Hamiltonian" 
in terms of the form factors in infinite volume becomes 



A{Lo) 


= Lo 


B{L,) 


= Lo 


C{Lo) 


= Lo 



rricLo rric 

2/cii 



:{A,m^{0)\A,{p)A,{-p)) 



m,Lo)3/2^ '' '^ ' '' ^^"^ ^^ '" ^^^ 



c^O 



^ :{AMM-p)\^mMp)M-p)) = %^ + ^ (3.16) 



where 



irricLof mlLo rric 



/n = F*(m,0) 



A{Lq) and B{Lq) does agree with the DGM result (3.13), but C{Lq) contains a new term 
coming from the disconnected part of the four-particle matrix element. The origin of this 
contribution is that in contrast to the bulk setting discussed in [1], the finite volume Hamil- 
tonian does not include any counter terms for the mass shifts (for details see Appendix 
C.2). The presence of the disconnected term is tested using TCSA numerical data in 
subsection 5.1.1. 

However, we already know that even (3.16) is inconsistent with the proper finite size 
description, which gives (3.11,3.12) to leading order in A. The underlying reason is that 
the density of two-particle states is not a product of the densities of one-particle states: 
the presence of the background phase shift modifies it by terms of order L~^ already at 
order A°. As previously, we continue to neglect residual finite size corrections (i.e. those 
that decrease exponentially with the volume), but all power-like corrections must be taken 
into account (at least to the necessary order in A) in order to maintain consistency. 

To determine the correct normalization of the states \Ai{pi)Ai{p2)) let us introduce the 
total momentum P = pi -\- P2 and the relative momentum p = {pi — P2)/2. The Jacobian 
of this change of variables is 1 and so (in infinite volume) 



{AM)MP2)\Mpi)MP2)) = {2^f5{P - P')5{p - p')yjpl + ml^pl + mf 

In finite volume, the density of states in the relative momentum variable can be obtained 
from the quantization condition (written here for the case of zero momentum ie. pi = 
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-P2 = P) 



pL + 5{E{p)) 
E{p) 



2n7r 



2\/ P + "^1 



(3.17) 



The density of states in the total momentum variable is unaffected, since the phase shift 
drops from the quantization of the total momentum, so it is only the density in the relative 
momentum of the two particles which must be modified. Let us examine the equation 

pL + 6{E{p)) = 2mT 

Taking the derivative of this equation with respect to p, first at fixed L and then at fixed 

n: 



L + 
L + 



dSjEjp)) 

dp 
d5{E{p)) 



271 — 

dp 

dL 

-p 



dp dp 

We obtain that the density of states at can be expressed as 

dn 1 dL(j)) 

dp 271 dp 



where the derivative 



is exactly the inverse of 



dL{p) _ dL 
dp dp 



dp{L) 
dL 
that was encountered in (3.10,3.11). 

Therefore the appropriate normalization of the zero-momentum two-particle state finite 
volume two-particle state is 



which gives 
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(3.18) 
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with p{L) given by (3.17). This is now consistent with (3.11) and the result of the Breit- 
Wigner analysis (3.12). The numerical TCSA results of section 5 do indeed confirm the 
inclusion of the "improvement" factors related to the density of states in finite volume. 
Finally, we wish to mention that similar "improvement" factors have already been proposed 
by Lellouch and Liischer for the measurement of weak decay matrix element of hadrons in 
[34]. 

4 Numerical methods 

4.1 Truncated conformal space approach 

In the framework of perturbed conformal field theory on a cylinder with spatial circumfer- 
ence L, the action (2.1) leads to the Hamiltonian 



H = H^ + fi dx^{t,x) + \j dx'^{t,x) 
Jo Jo 



where if* is the conformal Hamiltonian on the cylinder. Using Euclidean time r = —it 
and mapping the cylinder to the plane with 

27i{t — ix)\ _ f27i{T + ix 



z = exp ( J , z = exp 



we obtain 



27r I / - ex /L\^-'^* /-^ . . . , /^xi-2A* rL 



if = - I (^Lo + Lo - -j + /. ^-J J^ dxHz, -z) + A ^-J j^ «(., .-) 

(4.1) 
Due to translational invariance, the conformal Hilbert space Ti can be split into sectors 
characterized by the eigenvalues of the total spatial momentum 

Truncating these spaces by imposing a cut in the conformal energy, the truncated conformal 
space corresponding to a given truncation reads 

^Tcs(n,ecut) = ||^) e ^ I (Lq-Zo) 1^) =n|^), (^o + -^o- Y^j IV^) = e|^) : e < Ccutj 

The essence of the truncated conformal space approach as introduced by Yurov and 
Zamolodchikov [4] consists in realizing that on the space ?iTcs(^, Ccut) the Hamiltonian 
(4.1) becomes a finite matrix, which can be diagonalized numerically to get the finite 
volume spectrum of the model. 
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In the case of the Ising model (2.6), the Hilbert space can be written as 



' 16 ' 2 



where Vh is the irreducible Virasoro representation of highest weight h made from the 
Verma module 

\4 = {L-n, ■ ■ ■ L^rM : Lo\h) = h\h)} (4.2) 

by factoring out the singular vectors. Using conformal Ward identities^, the following 
matrices can be computed 

G^, = m 

(B^),, = {t\ail,l)\j) 
(5,)., = (^|e(l,l)|j) 

where \i) denotes a basis of ?iTcs(^? ^cut)- It is very convenient to use a basis which 
is directly related to the Verma module basis (4.2) which is, however, not orthonormal. 
Instead of performing a Gram-Schmidt orthogonalization procedure this problem can be 
remedied by considering the Hamiltonian 

which is isospectral to (4.1). We can introduce a dimensionless Hamiltonian measuring 
energy and length in units given by the lowest particle mass mi at the integrable point 

r = 

27r r /^^/^ 

hTcsA = —lHo + K,,—-—G-'B, + tKl/'''lG-'B, 

I = rriiL 

where we use the parameters and notations introduced in section 2.2. 

The application of TCSA to perturbations of c = 1 free boson CFT was developed 
originally in [7], and its use in double sine-Gordon theory is described in detail in [23, 27], 
therefore we give only a brief sketch to fix our conventions. For the case of the double 
sine-Gordon model (2.10), the relevant Hilbert space'' is 

H = 0^n (4.3) 

neZ 



^The algorithm we used here was developed originally for the numerical work in [6], but an explicit 
description of (the supersymmetric extension of) the matrix element determination can be found in [35]. 
The only published full TCSA algorithm [5] is not suitable for the present computation because it can 
only be appHed for low truncation levels, and is not automated enough to be scalable. 

''We consider only states with zero topological number, since the breather states we are interested in 
can be found in this sector. 
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where 

Tn = I a^ki---a-k„a_i^ ...a_ijn) : \n) = exp f n-v?(0,0) J |0) 

where the modes of the field are defined by the expansion 

ip{z, z) = (po 7=n log zz + -y= y] - (a„z"" + a„z~") 



Att V47r ^ n 



The matrices 

{Ho)ij = (a^ + A, 

(2|exp(nfv9(l,l))|j) 



T2) ^-■ 



{Vr 



njij 



vWMj) 



can be calculated explicitely in a closed form using the algebra of the free field modes. The 
dimensionless Hamiltonian then reads 

{ 2 4+3g 

(27r)T+f 2 (27r)^ 2 



where M is the soliton mass at the integrable point A = 0. 

In both models we only consider the sector with zero total momentum n = which is 
enough to obtain the necessary numerical data. The truncation levels in both cases were 
such that the Hilbert space dimension ranged from a few hundred up to 3-4000 states. 
This limitation was put mainly by machine time used for the numerical diagonalization, 
but in both cases the programs are fully scalable to any truncation levels, given enough 
computing resources (in terms of memory and execution time). For the Ising models 
we used Ccut = 19, ... , 27, while in the double sine-Gordon model the typical range was 
Ccut = 9, . . . , 15 for lower values of ^ and Ccut = 11, . . . , 17 for higher values (in one case we 
used Ccut = 18). In the case of double sine-Gordon theory, due to the choice S = —7i/2 it is 
possible to project the space of states onto even and odd sectors under the Z2 symmetry 

of the action (2.10) given by 

27r 

as described in [27]. This reduces the dimension of the Hilbert space with a factor of 
approximately 2 which makes it possible to use higher values for the truncation level. We 
used data from the even sector to extract the results as the choice of sector does not make 
much difference [27]. 
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4.2 Testing the numerics and the validity of form factor pertur- 
bation theory 

We tested the numerics by making several comparisons. First we took the integrable 
point, where we checked whether we get suitable agreement with the exact predictions of 
factorized scattering theory (masses, vacuum energy density and two-particle phase shifts). 
Next, in the case of the Ising model we checked the form factor perturbation theory 
predictions for the mass shifts (2.3) and the S matrix correction (2.4) using the form 
factors calculated by Delfino et al. in [2] (available at [32]) and formulae (2.3,2.4). This 
allowed us to determine the range of the coupling t in which perturbation theory holds 
with sufficient accuracy so that we can expect reasonable agreement for the perturbatively 
determined decay widths (2.9) as well. The S'-matrix comparison must be done carefully 
since formula (2.4) is valid at a constant center of mass energy. Therefore every formula 
must be converted from rapidity variables into the energy variable as follows. The phase 
shift can be extracted using 




^measured (£^) = -L^j {—j - mi(A)2 + 2mT 

where mi (A) is the mass of Ai including the correction to first order in A. This must be 
compared to 

<5thcorctical(^) = 5o{E) + 5i{E) (4.4) 

where the unperturbed phase-shift is 

/ E 

SJE) = -i log Su{^ = 2 acosh 

\ 2mi 

with mi denoting the mass at A = since this gives the energy dependence of the un- 
perturbed phase-shift, and to use (2.4) it is the energy that must be kept constant. The 
correction term can be written 



Si{E) 



where we can use again simply 



'd = 2 acosh 



^11 (^) 
E 



2m,i 
to first order in A. 

The agreement between the mass corrections and FFPT predictions for particles Ai 
and A2 is illustrated in figure 4.1, while the ^-matrix test is shown in figure 4.2. Both 
these tests indicated that |t| < 0.005 is a suitable choice for the range of the coupling. 

For the double sine-Gordon model at 6 = —7i/2 it is not possible to perform a similar 
test, because the mass and S matrix corrections vanish at first order in t and FFPT is not 
formulated to second order yet. However, the program was already thoroughly tested in 
previous numerical studies [23, 27]. We simply chose a range of the coupling in which the 
measured quantities were sufficiently close to linear dependence on t. 
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The mass of particle A1 



The mass of particle A2 
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Figure 4.1: Corrections to the masses of Ai and A2. TCSA data are indicated with crosses 
(the numerical uncertainties are too small to be displayed), while the lines give the FFPT 
predictions. 

4.3 Numerical methods for evaluation of the decay width 

We used three methods to extract the decay widths: 

1. The Breit-Wigner analysis (see section 3.1.2) gives the decay rate as 



mm5i{E) - max52(^) = 4:^/-pT 

where Si and ^2 are the phase shifts extracted from the two levels using (3.5) in 
the vicinity of the crossing. Because we aim to extract a/T to first order in the 
nonintegrable coupling t, rUa and rub in (3.5) must be substituted with the particle 
mass rrii to first order in t to keep consistency to this order. Furthermore 



P 



He 



<0 



E=Eo 



can be calculated using the t = exact S'-matrix for the background phase shift and 
Eq = rric the t = particle mass since any t 7^ corrections can be consistently 
neglected when comparing to the lowest order prediction (2.5) for the decay width 
r, from which we can extract the form factor value fdi- 

2. The "mini-Hamiltonian" method (section 3.1.1) when the matrix element B is ex- 
tracted from the minimum energy split (3.4) between the two levels 



5E{L, 



2\B\\ 



In a "naive" application of this method, we can then use the naive "mini-Hamiltonian" 
relations (3.16) (or, equally well, the DGM relations (3.13)) to link this directly to 
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Figure 4.2: The AiAi phase shift measured at t = 0, 0.005, 0.01 and compared to the 
prediction (4.4) as a function of the center of mass energy E (measured in units of the 
unperturbed particle mass mi). The effect of the resonance appears on the plot as a 
deviation from the first order FFPT prediction (2.4) in the range around Ec ^ rric ^ 
2.405mi, where a "shoulder" appears which can be observed quite clearly in the t = 0.01 
data. It is a fine resolution scan of this "shoulder" which is used to extract the decay width 
when applying the Breit-Wigner method. 
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the three-particle form factor: 

-BDGM(-^mm) 
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3. The improved "mini-Hamiltonian" method using the relations (3.18) derived in section 
3.2 ^ 

2 /cll 



B(W ' ^*<^) 



p dL 
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Here we replaced Lq by L^in which is valid to leading order in t. In addition, the 
state density correction factor can be computed from 

pL + 6t=o{E{p)) = 2n7r 

where 6t=o is the phase shift at the integrable point t = 0, again to leading order in 

t. 

In each of these approaches, the numerical setting is very convenient since the task is to 
find extrema of definitely convex or concave functions (the level splitting as a function of 
the volume or the phase shift as a function of energy) . Therefore it is not even necessary 
to find the level crossing volume Lq very precisely: any simple algorithm converges very 
fast even if the initial range given to start the search for extremum is quite wide. The 
background phase shift is given by (2.8) for the Ising and by (A. 5) for the double sine- 
Gordon model. Both give phase shifts monotonically decreasing with energy (in fact any 
two-particle scattering phases built as products of blocks {x} with < a; < 1 automatically 
possess this property). 

It is expected that the Breit-Wigner method should match the improved mini-Hamiltonian 
method as long as the resonance is narrow enough for the arguments of section 3.1.3 to 
apply. 

4.4 Systematic errors 

Contrary to the situation e.g. in lattice field theory, there are no significant sources of 
statistical errors. The TCSA Hamiltonian matrices are diagonalized numerically with 
double (16 digits) precision, and any error in calculating their spectra can be considered 
minuscule. There are, however, systematic errors in the extracted decay widths because (1) 
the TCSA Hamiltonian is not the exact one of the finite size system and (2) the analysis 
in Section 3 neglected residual finite size corrections (i.e. those that decay exponentially 
with the volume). 

The first source goes by the name "truncation errors", which originate from the fact 
that TCSA neglects an infinite tower of states lying above the truncation level, and grow 
with the volume. In the region close to a line crossing they manifest themselves most 
prominently in the fact that even for t = the degeneracy is generally not exact. This can 
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be modeled in the "mini-Hamiltonian" picture by adding a correction matrix representing 
the truncation errors to (3.43.2) 

As a result, the minimal splitting changes to 

6E{L^i,,) = 2\b + B\\ =2\B{X-Xo)\ (4.6) 

but B can still be read out from the slope of the dependence of SE{Ljain) on A and the 
form factor amplitude /di can then be determined from 

TJ 2/cii 



or 

B 



2/cii 



1 dp{L) 



p dL 



3/2 

L=Lo rric 



depending on whether we use the naive or improved "mini-Hamiltonian" method, respec- 
tively. 

This effect is demonstrated in figure 4.3, where it is obvious that the residual splitting 
decreases with increasing truncation level and gives an idea of how precise the effective 
description (4.5) really is. The effect of the other parameter a is to shift the value of Lmin 
while ^£"0 is the truncation error in determining the resonant energy Eq. 

It is easy to see that the same modification must be performed when applying the 
Breit-Wigner method, and relation (3.6) must be altered as 

mm6i{E)-m&x62{E) =4|5'(A-Ao)| 

from which the form factor amplitude /di for the process Ac ^ Aa + Aa can be read off 
using 

fell 



B' 




-^"^X/Mfe)^-1 
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In the sine-Gordon case, numerics always yields Aq = i.e. there is no residual splitting 
between the numerically determined levels at the integrable point. The reason is that the 
sine-Gordon spectrum is unaffected by A ^ —A since this can be implemented by shifting 
ip -^ ip + 271 /P in the action (2.10) and this symmetry is respected by the truncation 
procedure, so the numerical spectrum is even in A (up to very small errors in numerical 
matrix diagonalization). 

For a fixed value of Ccut and a choice of level crossing, we can then determine a value 
for fell- Then the next task is to extrapolate the results to infinite truncation level and 
then eliminate finite size effects (extrapolate to infinite volume). 
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Figure 4.3: The Breit-Wigner phase shift sphtting m.m6i{E) — max 62(E) as a function of 
the dimensionless coupling t in the Ising case for Ccut values of 11 and 29. 

Since truncation errors grow with increasing volume while finite size errors become 
larger with decreasing volume, some compromise must be struck in between. Masses are 
normally measured by finding the region where the measured gap is closest to being con- 
stant (has a minimum slope) and taking the value of the gap there. Two-particle phase 
shifts are measured over an extended volume range, and so they are affected by truncation 
errors at small, and by finite size errors at large values of the center of mass energy, and 
therefore it is very hard to perform any sort of optimization (except for a choice of the 
two-particle level used to extract the phase shift). Decay widths are measured in small 
regions around level crossings, and we perform the measurement on the first few of them, 
for several different values of Ccut • 

Truncation errors can be improved by fitting the truncation level dependence by some 
function for extrapolation to Ccut = 00. Detailed examination of numerical data shows that 
they can be fit in most cases by a function of the sort 



f{e,^t,L) =f{L) + a{L)e 



X 

cut 



(4.7) 



where f{L) is the extrapolated value of the measured physical quantity / at volume L, 
and the second terms describes the truncation effects with a and the exponent a; > 
to be determined from the fit. Such an extrapolation was used in [36] where it proved 
very useful in the determination of vacuum expectation values. Experimenting with this 
technique in cases where a theoretical prediction for the measured quantity is available 
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Figure 4.4: Truncation level extrapolation for s^n in the double sine-Gordon model (i? is 
related to i via 2R^ = 1 + i'^) 

shows that this really helps, and in many cases the truncation errors are improved by an 
order of magnitude. Based on the numerical results the exponent x can be supposed to be 
independent of the coupling t and depends only very mildly on L (in the sine-Gordon case, 
however, a dependence on the parameter ^ must be taken into account). For an illustration 
see figure 4.4 which shows that separate extrapolation is necessary for odd and even values 
of the truncation level. The reason is that only when increasing the truncation level by 
a step of 2 we get new vectors in each (Verma or Fock) module: since the momentum is 
zero, we must increase the descendent quantum numbers on both left and right by 1, and 
thus it is only consistent to group together data pertaining to Ccut that differ in steps of 2. 
In the sine-Gordon case the situation is further complicated by the fact that the number 
of Fock modules JF„ contributing states in (4.3) also grows with the truncation level, but 
this effect turns out to be numerically irrelevant. 

Unfortunately, there is no theoretical method providing a reliable estimate of trunca- 
tion errors. However, in the case of the double sine-Gordon model the two Breit-Wigner 
extrapolation curves approach the average value from the two sides (as in figure 4.4 (b)) 
and therefore residual truncation errors can be estimated by the difference between the two 
extrapolated values. We adopt the same procedure for the "mini-Hamiltonian" methods 
as well, although in that case the situation is less clear-cut because the two sets of data 
approach the extrapolated limit from the same side. In the case of the Ising model, the 
even/odd extrapolations approach the final value from the same side for both methods, so 
the difference between is less reliable as an estimate for the truncation errors. 

The other source of systematic errors is residual finite size corrections originating from 
vacuum polarization and the finite range of interactions. According to Liischer's result [30, 
31] these are suppressed by a factor of the form exp(— ML) where M is some characteristic 
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mass scale (typically the mass of the lightest particle). Such effects are e.g. the volume 
dependence of the particle masses and vacuum (Casimir) energy. These were neglected in 
all the considerations in Section 3, where the only finite size effects taken into account are 
the ones that decay as a power (generally L'"^) of the volume. 

Having taken care of truncation errors, residual finite size corrections can in principle 
be suppressed by fitting an appropriate extrapolation function. Suppose that the decay 
width can be calculated from several crossings which take place in very small regions around 
different values of the volume Li, . . . , L„ and that truncation errors in these determinations 
can be kept suitably small. Using our general understanding of finite size effects, we can 
think of the decay form factor /di determined at a crossing at Lj as the value of a function 
fcii{L) at L = Li with the property that 

/cii(^) = /cii + 0(e-^^^) 

where /di is the value predicted by the infinite volume quantum field theory. Unfortunately 
our numerical results (even after extrapolation in the truncation level) do not have sufficient 
precision to fit the exponential term in a reliable way. Therefore we will not perform this 
extrapolation in the sequel, and restrict ourselves to quoting the values of fdi as measured 
at each level crossing separately in the Ising case; in the sine-Gordon case truncation errors 
allowed us to use only the first level crossing, i.e. the one that occurs at the smallest value 
of the volume. 

5 Results 

5.1 Ising model 

5.1.1 Testing the mini-Hamiltonian coefficient C 

It is interesting to check whether the improved "mini-Hamiltonian" description (3.18) is 
consistent with the numerical data. The correct formula for B can be tested by comparing 
the results for the decay width extracted using the Breit-Wigner method, and will be 
preformed in subsection 5.1.2 below. Here we perform a quick test for the coefficient C. 
Let us take the sum of the two energy levels Ei{L) + E2 (L) at the volume where the level 
crossing occurs at A = 0. Then using the results of subsection (3.1.3) this must be equal 
to 

Ei(Lo) + E2iLo) = 2m, + A(A(Lo) + C(Lo)) 

when taken at L = Lq where Lq is the volume where the levels cross. Using TCSA data 
with truncation level Ccut = 27, we measured this energy at five different values of the 
integrability breaking coupling t = —0.003, —0.001, 0, 0.001, 0.003, at the the first three 
level crossing. The results were then fitted by a linear function, parameterized as 

En{Lo) + EciU) 

= a + bt 

mi 
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Numerics (TCSA) 


"mini-Hamiltonian" prediction 


Lo 


a 


h 


h (DGM) 


b (naive) 


b (improved) 


18.152 


4.52877 


-3.598 


-0.509 


-3.580 


-3.974 


24.900 


4.80039 


-4.348 


-0.947 


-4.016 


-4.364 


34.184 


4.81197 


-4.586 


-1.263 


-4.334 


-4.616 



Table 5.1: Testing the effect of the disconnected term and the finite volume state density 
factor on the "mini-Hamiltonian" description. We compare the TCSA data to three al- 
ternatives: the DGM result (3.13), the naive "mini-Hamiltonian" (3.16) and the improved 
version (3.18). Due to residual finite size effects, the comparison between the naive and 
the improved "mini-Hamiltonian" at the first level crossing cannot be taken seriously. 

and are summarized in table 5.1. The theoretical value of a is 

2m4 



4.80973 . . . 



m\ 



It is obvious that to get agreement with the TCSA results, both the disconnected 
term and the correct finite volume normalization of two-particle states must be taken into 
account. Another important point is that residual finite size effects are largest at the first 
level crossing, where a = 4.52877 is quite far from the theoretical value, therefore one 
cannot distinguish between the naive and the improved "mini-Hamiltonian" methods. On 
the other hand, the results at the other two level crossings show that both the disconnected 
terms and the improvement is needed in order to get appropriate agreement with the 
numerical data. 

5.1.2 Decay widths 

We measured decay width for the process A4 ^' Ai + Ai at three level crossings, while 
Aq -^ Ai + Ax was measured at one level crossing. We extracted the form factor matrix 
elements /411 and /511 using all the three methods ("naive" mini-Hamiltonian - equivalent 
with the DGM one for the case of the decay width-, improved mini-Hamiltonian and Breit- 
Wigner method) and the end results (after truncation level extrapolation) are presented 
in tables 5.2 and 5.3. 

From table 5.2 we can see that indeed it is the improved "mini-Hamiltonian" method 
(rather than the naive one) which is consistent with the Breit-Wigner approach. It is 
also apparent that the extracted values do change with the volume, and that the naive 
"mini-Hamiltonian" method significantly overshoots the theoretical prediction, while both 
the improved "mini-Hamiltonian" and Breit-Wigner methods are in reasonable agreement 
with it (within 1% at the 3rd level crossing). 

In the case of the matrix element /511 truncation errors are expected to be higher 
because the levels lie much higher in the spectrum. This observation is borne out by the 
data presented in table 5.3, where we gave the data at the two largest values of Ccut instead 
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n 


miLo 


naive mH (o/e) 


improved mH (o/e) 


Breit-Wigner (o/e) 


1 


18.152 


37.658/37.662 


33.422/33.426 


33.255/33.336 


2 


24.900 


38.871/38.939 


35.736/35.799 


34.574/34.887 


3 


34.184 


39.099/* 


36.829/* 


36.318/* 



Table 5.2: The measured values of /411 where n is the label of the level crossing, Lq is the 
place of the level crossing for t = and the last three columns give the result extracted 
using the method indicated at the top. The two numbers given are the extrapolated values 
in the truncation level for Ccut =odd/even (*s indicate the cases where the data could not 
be reliably extrapolated). The theoretical prediction is /411 = 36.730. 



n miLo naive mH (27/28) improved mH (27/28) Breit-Wigner (27/28) 


1 23.206 21.011/21.120 20.244/20.349 18.551/18.612 



Table 5.3: The measured values of /511 at the first level crossing. Lq is the place of the level 
crossing for t = and the last three columns give the result extracted using the method 
indicated at the top. The theoretical prediction is /511 = 19.163. Data at the second 
and higher crossings already have too large truncation errors. Some of the data could not 
be extrapolated in the truncation level; therefore we only quote the values measured at 
Ccut = 27 and 28. 

of extrapolating them in the truncation level, because in most cases the extrapolation 
function (4.7) could not be fitted in a reliable way. The variation of the primary quantities 
A6 used in the Breit-Wigner method and the minimum energy split used in both mini- 
Hamiltonian methods as a function of t also shows much larger t^ effects (but decreasing t 
is not possible beyond a certain limit because of the presence of the residual line splitting 
4.6). Therefore the matching between the Breit-Wigner and the improved "mini-Hamilton" 
methods is less precise, but the latter is still somewhat closer to the Breit-Wigner result 
than the naive version. Despite the much higher errors the two best estimates 20.349 (from 
the improved "mini-Hamiltonian" method) and 18.612 (from the Breit-Wigner method) are 
still within a few percent of the theoretical value 19.163. 

5.2 Double sine-Gordon model 

For convenience, we parameterize the sine-Gordon coupling as follows 

/4^ . 1 



P 



R 



e 



2i?2 



In table 5.4 we compare the measured values of S311 (see eq. (2.13)) for the three 
methods with the predictions, using the first level crossing (for a plot, see figure 5.1). For 
the "mini-Hamiltonian" methods, the values in the table are the averages of the two sepa- 
rate extrapolations for Ccut =even/odd, while the errors shown are the estimated residual 
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Figure 5.1: The theoretical value of s^u plotted against ^, together with the numerical 
results obtained with the naive and the improved "mini-Hamiltonian" methods (Breit- 
Wigner data are omitted, since they are not significantly diff'erent from the latter). The 
error bars include only the residual truncation errors (the residual finite size errors are not 
estimated). 
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R 


nmH 


imH 


BW (e/o) 


FFPT 


ML, 


g-miLo 


1.44 


0.931 ±0.025 


0.821 ±0.021 


* /* 


0.7759 


16.469 


10-^ 


1.5 


1.328 ±0.038 


1.185 ±0.034 


* /* 

/ 


1.1694 


12.398 


2 ■ 10-^ 


1.6 


0.959 ±0.012 


0.883 ±0.011 


0.904/0.893 


0.9303 


11.588 


0.0002 


1.7 


0.630 ±0.004 


0.592 ±0.003 


0.584/0.587 


0.6383 


11.927 


0.0005 


1.9 


0.286 ±0.001 


0.274 ±0.001 


*/0.269 


0.2917 


13.615 


0.0011 


2.2 


0.1076 ±0.0002 


0.1046 ±0.0002 


0.110031/0.110031 


0.0999 


17.475 


0.0018 


2.5 


0.04867±2-10-5 


0.04767±2-10-^ 


*/0.049951 


0.0392 


22.309 


0.0023 


2.6 


0.03845 ±6 -10-5 


0.03773 ± 6 ■ 10-5 


0.041271/0.041283 


0.0295 


24.089 


0.0024 


2.7 


0.03075 ±5 -10-5 


0.03022±5-10-5 


0.033467/0.033403 


0.0224 


25.946 


0.0025 



Table 5.4: Measured values of S311 using the "mini-Hamiltonian" and Breit-Wigner meth- 
ods. The column label "nmH" corresponds to the naive, while "imH" to the improved 
"mini-Hamiltonian" method. The next column contain the results of the Breit-Wigner 
method extrapolated for even/odd truncation levels, respectively, which also give a rela- 
tively good estimate of residual truncation errors. FFPT labels the theoretical prediction, 
and the last two columns contain the volume corresponding to the level crossing at t = 
and the characteristic suppression factor of the leading finite size correction. The *'s la- 
bel the cases where the Breit-Wigner data could not be extrapolated meaningfully in the 
truncation level. 

truncation errors calculated dividing the difference of the two extrapolated values by 2. 
Because the truncation level extrapolation for the Breit-Wigner method is more problem- 
atic (at least for some values of i?), we show the results for even/odd truncation levels 
separately. Truncation errors increase with decreasing R, and indeed this is most promi- 
nently shown by the relative deviation of the measured and theoretical value at i? = 1.44 
where the Breit-Wigner method is not applicable due to the fact that albeit theoretically 
the condition of the monotonous decrease of the background phase shift 






<0 



is still satisfied, the numerical data fail to have this property. In addition, the value of 
Lq for R = 1.44 is also significantly larger than for R = 1.5, which increases truncation 
errors too. For larger R the deviation between TCSA numerics and FFPT seems to be 
controlled by the residual finite size effects. Higher level crossing show too large truncation 
effects, and are not useful for numerical analysis. Therefore there are not enough data to 
estimate the residual finite size effects, and to get some idea of their magnitude we quote 
the exponential suppression factor characterizing their decay with the volume. Note that 
this is not a very precise estimate since there must certainly be additional dependence on 
^ and on some power of Lq, the form of which is not known. However, it captures the 
tendency correctly: the residual finite size corrections increase with R, again making the 
measurements at higher values of R less precise. 
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It is also obvious that the improvement of the "mini-Hamiltonian" method is most 
important for smaller R and for R = 1.5 . . . 1.9 it is required to get appropriate agreement 
with the Breit-Wigner method. For larger values of i? the improvement is too small to draw 
any useful distinction between the naive and the improved methods. At R = 1.44 we see 
that despite the problems with the Breit-Wigner method, the improved "mini-Hamiltonian" 
method is still applicable and is in reasonable agreement with the theoretical prediction 
(relative deviation 6%), while the naive (DGM) version is again off the mark with a relative 
deviation of 20%. 

On the other hand, note that the deviations in the range R = 1.44 . . . 2.2 obtained at 
the first level crossing are comparable in magnitude to those obtained in the case of the 
Ising model for the process A4 -^ Ai + A\ at the second level crossing and so we consider 
the agreement with the theoretical predictions fully satisfactory. 

6 Conclusions and outlook 

In this paper we investigated the signatures of resonances in finite volume quantum field 
theory. The central result of the paper is the development of two methods: (a) the Breit- 
Wigner method utilizing the standard parameterization of the resonant contribution to the 
two-particle phase shift, and (b) the improved "mini-Hamiltonian" method based on the 
explicit Hamiltonian description of quantum field theory, to extract (partial) decay widths 
of particles given the finite volume spectrum. 

We have shown the equivalence between these methods to first order in the resonance 
width (i.e. for narrow resonances). As a further result of this investigation, we now 
have a finite size description of the resolution of level crossings which brings Liischer's 
description of finite size effects, the Breit-Wigner parameterization of resonances, the 
"mini-Hamiltonian" method and form factor perturbation theory together in a consistent 
framework. 

Both methods are extremely simple to implement numerically since all what is needed 
is to find a minimum/maximum of some convex/concave function (the level splitting as a 
function of the volume or the phase shifts extracted from the two levels as a function of 
energy). However, the Breit-Wigner method does not seem to be as robust numerically 
under extrapolation in the truncation level (i.e. the UV cut-off) as the improved "mini- 
Hamiltonian" method, so the latter appears to be preferable. 

Using the machinery available in the realm of 1+1 dimensional quantum field theo- 
ries, namely form factor perturbation theory and truncated conformal space approach we 
demonstrated the efficiency of these methods, by showing agreement between the numer- 
ically extracted decay matrix elements and the theoretical ones in many cases up to a 
precision of 1 to 7% (as shown in figure 5.1). The most important factor in achieving this 
level of precision this is that both methods measure directly effects proportional to VT 
(where F is the resonance width) in contrast to other methods proposed in the literature 
(such as the "plateau slope" method in [8]) which attempt to detect effects of order F. 
Therefore, our methods are much more sensitive when the resonance is narrow. Another 
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important advantage of methods (a) and (b) is that they only need data from a small 
vicinity of the level crossings resolved by the finite width of the resonance, while (as dis- 
cussed at the end of subsection 3.1.2) fitting a Breit-Wigner shape directly to the extracted 
phase shift is not as reliable because it uses data from an extended range of the volume 
(containing two subsequent level crossings), and the residual finite size effects introduce 
significant distortion in the shape of the resonance curve. 

It is clear that since the necessary description of finite size effects exists in arbitrary 
dimensions [31, 30, 8], the development of the two methods in Section 3 can be generalized 
to any number of space-time dimensions in a straightforward way and the results of the 
present work can be extended to more interesting theories formulated e.g. on the lattice. 
We remark that even the form factor approach to the calculation of the decay width has its 
3 + 1 dimensional counterpart in the description of the weak decays of hadrons, where the 
appropriate form factor is the hadronic matrix element of the weak current. The crucial 
issue is whether the necessary precision in the measurement of finite size effects can be 
attained. Therefore the most important open problem is to extend these results to theories 
in physical (3 + 1) space-time dimensions, and in particular to phenomenologically relevant 
models. 

Acknowledgments 

G.T. would like to thank Z. Bajnok for useful discussions and for drawing attention to 
paper [2]. This work was partially supported by the EC network "EUCLID", contract 
number HPRN-CT-2002-00325, and Hungarian research funds OTKA T043582, K60040 
and TS044839. G.T. was also supported by a Bolyai Janos research scholarship. 

A The 311 form factor 

In this paper we only calculate the decay width for the simplest possible case B^ -^ 
Bi + Bi. In order to accomplish this, we need the form factor of the perturbing field 

^ = cos (1$ + 6) 

i^3n (^1,^2,^3) = {0\cos(J^^ + 5] |i?3(^i)i?i (^2)^1(^3)) 

In [37] Lukyanov obtained a closed expression for the form factors of exponential operators 
with any number of Bi-s in the asymptotic state. It can be written as an expectation value 

Fl, (^„, . . . , ^1) = (0| exp (m$) \B,{^^) . . . B^{A)) = GaiH^n) ■ ■ ■ HA)) 

where 



^ ' 2sin7rf V / 




A (4) = 2 cos I -— W 2 sm I -— exp 



2 y y V 2 / \ Jo 27r sin t 
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and the expectation value can be computed using Wick's theorem for the generalized free 
field uj (in a multiplicative form) with 



N N-l N 

j=l j=l k=j+l 



\OlhOlj 



where 



R{^) = Arexp<^8 



Af = exp<A 



(it sinh t sinh t,^ sinh t(l + ,^) 



t sinh^ 2t 

dt sinh t sinh tC, sinh t(l + ^) 
t sinh^ 2t 



i^ 



sinh^ t I 1 

TT 



(A.l) 



(the integral representation is only valid in the strip — 27r + n^ < '^m'd < — tt^) is the 
so-called minimal -Bi-Bi form factor satisfying 



R{-m + d) = Ri^-m-d) 
where »S'ii(^) is the Bi — Bi two-particle scattering amplitude. Furthermore 



(A.2) 



ga{P) = {e^n 



87r-/32 



2r 



X exp 



/3V2 
87r-/32 

°°rft 

T 



(A.3) 



-nh^ (if ^) 



-2t 



2 sinh ( ^t ) cosh f f 1 - g ) t ) sinh t 



47r 



is the exact vacuum expectation value of the exponential field [38]. 
Performing the Wick contractions explicitly we obtain 



(A(^i)...A(^„)) 



A 



2 sin vr^ 



En 



ajC ' ~ I X 



n—l n 



WWRi^k- ^j - {Ck- a,)i\ 



7l\°'kOl3 



j=l k=j 

where the sum goes over all the 2" possible "configurations" of the auxiliary variables 
aj = ±1. Using the identity 



R{d)R{i9 ± in) 



sinh -(9 



sinh {} ^ i sin nC^ 



(A.4) 



38 



we can write 



-ji-^akOj ^^ ^ ,ak — Oij siniT^ 



The end result is 

n—l n 



V 2sm7rey J-^A^Al^ 

The parity properties of the form factors are determined by 

p-a __ I 1 \n pa 

n \ J n 

As a result we obtain 

(0|cos(a$)|fii(^i)...fiiK)) = ^-^(0|exp(2a$)|5i(^i)...5i(^„)) 

(0|sin(a<l>)|5i(^i)...5i(^„)) = i^lzil!(o|exp(2a<l>)|5i(^i)...5i(^„)) 

Form factors of higher breather can be computed using the bootstrap fusion rules [9]. To 
obtain B^ one can fuse 2 _Bi-s into a B2 and then a Bi and a B2 into a B^. The fusion 
angles and three-particle couplings can be read off the Bi — Bi and Bi — B2 S matrices 
[29] 

sinh {} + isirnif . . x 

-^11 = —r-. ^ A.5 

smh V — t sm vr^ 

sinh {} + isin ^f sinh 'd + isin ^C 

sinh i!) — isin |^ sinh 1!) — isin ^^ 

If the particle c occurs as a bound state of a and b then the corresponding S matrix has a 
pole 

•| c 12 



^ - ^"afe 



from which we obtain (choosing the three-particle couplings real and positive) 



uli = TT^ , 7ii = a/2 tan n^ 



o Stt q / tanTT^tan^^ 

Therefore the 311 form factor can be obtained starting from a 5-particle 11111 form factor 
and fusing the first three Bi particles into a B3. This means that the decay process 

39 



B^ ^' Bi + Bi has a non-vanishing amplitude to leading order only if (5 7^ 0. Indeed, for 
6 = this process (and its generalization -B2n+i ^ -Bi + -Bi) is entirely forbidden by the 
C-parity of the breathers. 

Using the notation Xj = e^^ we get the following result for the 5-particle form factor 



F,\,,,{A,...,A) = QaX'Qtuu{A,...,A)l[ n Ri^k-^j 



(A.6) 



j=l k=j+l 



Ql,,,, {§,,...,§, 



, , , , 14 r 12 , ([a]^(fflf74-fX5) +f^5)(f^2fX3-4(75COs2 7rO ^ 

—i\a\ < \a\ — \a\ -\ > 



n n (^i + ^fc) 
j=i k=j+i 



where the elementary symmetric polynomials (Jk are defined by the generating function 

5 5 



i=i 



fc=0 



and 



sm 



TT^a 



sinTT,^ 

In the following we choose 5 = —7i/2 because it will be convenient later due to another 
Z2 symmetry which exists at this point (and which, on the other hand, forbids every 
process B2n -^ Bi + Bi to the lowest order) . Then the form factors of the perturbing field 



^ = sin 1$ can be written as 



P. 



/3/2 



Fl, (^1, . . . ,^„) = (0| sin {^<!>\ |i?i(^i) . . . i?i(^„)) = -tFl[t (^1, . . . ,^„) 

for n odd. 

First we fuse two -Bi's to obtain the form factor F21U. The relevant equation is 



^7ni^2ni (^1,^3,^4,^5) = ResF*iii U + ^ (iTv^ - e) ,^1 - ^ (m^ 



e),^3,^4,^5 



The only singularity comes from the factor i?(^2 — ^1) in (A.6). The residue can be 
calculated quite simply due to (A. 2) with the result 



Res R (^) = Res ^n (-^) R (-^) = i hfiVR («<) 

1?=— ivr^ 1?=— ivr^ 



and we get 



Fl,,{A,^3,KA) = ^;3/2AH\^(mOgfnilh?l + ^,^l-^,^3,^4,^5 X 



i< 



in^ 
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4 5 



n n R{^k-^j)i[R(^i-A 

j=3 k=j+l 1=3 

where Qum = -iQuui 



^"^^i?U-^, + ^"^ 



Now we fuse Bi and B2 into -B3 



^l!2Fi^^ (^1,^4,^5) = ResF^, (^^, + ^ 
Using the same tricks as above plus the identity 



2^< - 2 ' ^1 ~ ^^^ + 2 ' ^^' ^^ 



, tan ^ 

5ii (2^0 - 2 



tan 



TT^ 



3 \ 2 

7?i 



we obtain 

Qniii (^1, ^4, 171^ + ^5, ^5, -^vre + ^5, ) X 

R (^^5 - ^4) n [^ (^i - ^1 + ^^0 R i^j - ^1) ^ (^i - ^1 - ^ttO] 

j=4,5 

To obtain the decay rate we need 

/311 = |i^3n(^^>^c,-^c)| 

where the rapidity ^c of the outgoing Bi particles can be computed from 



mi sm ^ 



The result can be written as 



where 



/an = Qm>^'lnlt2Qnni (0 ^sn (0 



Qiiiii(0 = -gfnii(^c,-^c,^vr(l + 0,^vr,^7r(l-0) 
_ (l + 2cosO(l + 2cos< + 2cos27rO 

64 cos TT,^ cos^ ^ 



and 



7^3ll (0 



Riyinif R{2mi) R{-2d^) x 

R (^e - ^7r(l - 0) i? (^e - ^vr) i? (^e - ^7r(l + 0) X 
i? {-de - mil - 0) R (-^c - ^vr) i? (-^e - ^7r(l + 0) 
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We can make use of eq. (A. 4) to shift the arguments of the i?(2^c), -R(^7r^) and R{2in^) 
factors into the validity range of the integral representation (A.l). The result is 



7^311 (0 



sinh 2'dr 



cosTT^ R (^c - i7r{l - Of R (^c - iT^f R {^c - i7r(l + ^f 



sinh 2^e - i sin < 2 cos< + 1 R (i7r({ - 1))^ R (i7r(2^ - 1)) R (-2^^ 



tTT] 



For the dimensionless matrix element ssn we obtain 



S311 (0 



/; 



311 



M2+2« 



^(OA(0'2tan7re. 



'tan 



Sire. 



tanf 



■Qiim (0^311 (0 (A.7) 



where 



^(0 



2r(|) 



2+2^ 



exp 



°° dt 
2t 



sinh (^) 



cosh ( YT7 ) sinh t '- ~^ ^ 



e g-2t 



B Extracting A6 = min Si—insixS2 from the "mini-Hamiltonian 

We start from the expression (3.7) for the two levels 
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-^1,2 = mc + A\- a{L - L^nn) ± v «^(-^ - -^ 



\2 I D2\2 



52A2 



From this the phase shift functions 61^2 defined in (3.5) can be expressed as functions of L 



h,2{L) = -L^ 



EiAL) 



mt 



We demonstrate how to find the extremal value of Si{L). We solve the equation 

d6i{L) 



dL 







by substituting L = Lmin + k\ and expanding in A. The lowest non- vanishing order is A° 
and we find the following equation for n 



ml — Ami ~ LqiiIc ( a 



a^K 



VWT 



a^K^ 







Solving this equation and substituting the solution back into (5i, we find that to linear 
order in A 

a/— (m^ — Am\ — 2aLQmc) 



min(5i = -Lo^/(^^) ~mj + B\ 



2a 
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Performing a similar calculation for S2 we obtain 



max52 = -LqW f — ^j -mf-BX 



mc\2 \/-{jnl - 4"^i - '^OiLomc 



2a 
The difference is 

mmoi — maxo2 = -dA-^^ (B-l) 

We can go further by calculating the slope of the phase shift at A = 0. The energy of the 
two-particle level (the one which is not independent of L) is then 

E = rric - 2a{L - Lo) (B.2) 

Using 

E-' 



we get 



6,^o{E) = -L{E)\l[-] -ml 



d6x=o 



dE 



ml — Arnl — 2aLom,c 



E=mc ^a^ml - Amj 

which must be negative for the Breit-Wigner method to apply, so eq. (B.l) makes sense. 
To leading order in A Sx=q can be substituted with the background phase shift ^o and m,c 
with the resonance position Eq and so we obtain 



E=Eo 



min(5i — max(52 = 23X^1 y^ — 4mf —— 

C The "mini-Hamiltonian" coefficient C 

C.l Determining C from finite size corrections 

Let us calculate the shift of a two-particle level to first order in A. For zero total momentum, 
the relations (3.1) can be reformulated as 



L\l{^y -m,{Xy + 6{E{X),X) = 2nn 

To first order in A, we get the following equation for the energy shift SE: 

^ ^-m,5m, ^ d6o{E) ^^ ^ dS{E,X) ^^ ^ ^ 
' ^2 Z dE dX 



(f) -^ 



2 
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where E is the energy, mi is the particle mass and So{E) is the phase shift at A = 0, which 
satisfy 



^v(f) -^i + ^o{E) = 2n7r 



(C.l) 



and Srrii is the mass shift to first order in A. As a result 

Lmi6mi - 



(!) 



2 o d5(E,X=0) . 



5E 

¥ + \/(f)'-'"?^ 

Using (2.4) we can write at the level crossing L = Lq (recall that E{Lq) = rric) 
mr\^ ^2 d6{m„ A) ^ iQ-^Mmc)Jf!!h 

A=0 



— m 



d\ 



mi 



dSii{mc,\) 



dX 



A=0 



nil siiih i^i 



(cll) 



5n (-2^f ')) F*ii (tir, 2^f '^ + ^tt, 0, 2^f'') 



7* L()(cll) 
1111 



m^ sinh 2'i9-['^ 



m^ 



From (2.3) 



rriiSmi = AF*(i7r, 0) 



On the other hand, using E{L) = 2^Jp{L) + m\ (with the two particles having momentum 

±p(L)) we get from (C.l) 



LE 

— + 

and so we obtain 

1 dp{L) \ 



E 

2 / ^ dE 



,dSo{E) pEdL 



rrir 



L=Ln 



C 



p dL J 



4 dp L=Lo A ( 1 dp{L) \ 
\ P dL J 



u 



L=Lo 



mt 



rrir 



C.2 Disconnected parts of the four-particle matrix element 

We need to determine the matrix element {Ai{p)Ai{—p)\'^{0)\Ai{p)Ai{—p))L=Lg. Return- 
ing to rapidity variables and using the crossing properties of form factors we can write (in 
infinite volume) 

i„(Ai(^3)A(^4)|^(0)|A(^i)A(^2))i„ = F*ii(^3 + ^vr,^4 + ^7r,^i,^2) + 

27r5(^i-^3)i^ii(^4 + ^vr,^2) + 
27r5(^2-^4)i^ii(^3 + ^7r,^i) + 
27r5(^i- ^3)5(^2- ^4) (0|^(0)|0) 
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The last term can be dropped because it is related to the vacuum energy shift, but we 
normalized the finite volume energy levels by subtracting the ground state. In [1], the 
other two disconnected pieces were canceled by mass shift counter terms; however, the 
finite volume Hamiltonian does not contain such terms. The S functions can be written as 

2tt6 {'di — '1)2) = 2'n'6{pi — p2)iTii coshi^i 

and in finite volume (Lq) 

2nS (^1 - ^2)1^1=^2 = Louii cosh^i 

We need the matrix element at the special rapidity values ^1 = —^2 = ^3 = —^4 = ^1 
(recall that rric = 2mi cosh.i9]^ ) which gives 



J 

nicLoF^^ {in, 0) 
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